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ABSTRACT 

Observations of distant galaxies play a key role in improving our understanding of the 
Epoch of Reionization (EoR). The observed Lya emission line strength - quantified 
by its restframe equivalent width (EW) - provides a valuable diagnostic of stellar 
populations and dust in galaxies during and after the EoR. In this paper we quantify 
the effects of star formation stochasticity on the predicted Lya EW in dwarf galaxies, 
using the publicly available code SLUG (used to 'Stochastically Light Up Galaxies'). 
We compute the number of hydrogen ionizing photons, as well as flux in the Far UV 
for a set of models with star formation rates (SFR) in the range 10~ 3 -1 M yr _1 . 
From these fluxes we compute the luminosity, L a , and the EW of the Lya line. We 
find that stochasticity alone induces a broad distribution in L a and EW at a fixed 
SFR, and that the widths of these distributions decrease with increasing SFR. We 
parameterize the EW probability density function (PDF) as an SFR-dependent double 
power law. We find that it is possible to have EW as low as ~EWo/4 and as high 
as ~ 3xEWo, where EWo denotes the expected EW in the absence of stochasticity. 
We argue that stochasticity may therefore be important when linking drop-out and 
narrow-band selected galaxies, when identifying population III galaxies, and that it 
may help to explain the large EW (EW >100 — 200 A) observed for a fraction of Lya 
selected galaxies. Finally, we show that stochasticity can also affect the inferred escape 
fraction of ionizing photons from galaxies. In particular, we argue that stochasticity 
may simultaneously explain the observed anomalous ratios of the Lyman continuum 
flux density to the (non-ionizing) UV continuum density in so-called Lyman-Bump 
galaxies at z = 3.1, as well as the absence of such objects among a sample of z = 1.3 
drop-out galaxies. 
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1 INTRODUCTION 

The Epoch of Reionization (EoR) represents a milestone in 
the evolution of our Universe, during which hydrogen gas 
in most of its volume changed from fully neutral to fully 
ionized. Observations of the cosmic microwave background 
support that reionization started at redshift z ~ 20 while 
quasars observations demonstrate that the process was com- 
pleted at z ~ 5 — 6 (e.g. Mesingcr 20l"o| |Pritchard et al.| 
20101. One key to new insights in this area are observations 



of high redshift galaxies, as they will allow for the construc- 
tion of a complete physical picture of the reionization pro- 
cess. 

* IAU Gruber Fellow, Email: forero@berkeley.edu 



Observations of distant Lyman Break Galaxies (LBGs) 
and Lyman a Emitters (LAEs) in the redshift range 4 < 
2 < 10 have provided the most prominent windows on the 



star formation history at these early epochs (see e.g. Robert- 



son et al. 2010 and references therein). LBGs are galaxies 



that have been selected using broad-band surveys, which 
are sensitive to the 'Lyman break' in the spectra of star 
forming galaxies. This break is a consequence of typical ef- 
fective temperatures of O & B stars, and of interstellar and 
intergalactic absorption of flux blueward of the Lya wave- 
length (Steidel et al. 19961. Surveys for LBGs (or 'drop- 



out' galaxies) are most sensitive to UV rest-frame contin- 
uum emission of galaxies. LAEs are galaxies that have been 
selected from narrow-band surveys which are sensitive to 
high-redshift galaxies that have strong Lya line emission 
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(e.g. |Rhoads et al.|2 000). The Lya line originates in HII re- 
gions that surrounded O & B stars. As a result, the Lya line 
flux (and strength) in LAEs are more directly connected to 
the ionizing photon budget (e.g. |Schaer cr 2003). 

The 'strength' of the Lyman a emission line is quanti- 
fied by its rest-frame equivalent width (EW), which corre- 
sponds to the ratio between the line intensity and the contin- 
uum. The Lya EW is an important physical quantity associ- 
ated with galaxies during the EoR for (at least) three reasons 
: (i) narrowband surveys for LAEs are typically sensitive to 
galaxies for which the Lya EW exceeds some threshol d value 

65 A 



sec e.g. 



Ouchi 



(which can range between EW~ 20 
et al.|[2008[ ). The Lya EW is clearly an important physical 
quantity that links the LBG and LAE populations; (m) the 
'first' galaxies that formed from pristine gas are expected 
to be characterised by unusually large EW (well in excess 
of EW~ 200 A), which may be one of their most promi- 
nent observational signatures ( Schaerer|2003 Johnson et al. 
2009 Raiter et al.||2010 i; (Hi) recent observations of LAEs 
& LBGs in the redshift range z = 3 — 7 have revealed an 
apparent non-monotonic evolution in the EW distribution, 
with significantly less Lya being observed from galaxies be- 



yond redshift z > 6 (Stark et al.|2010 


2011 


Pentericci et al. 


2011 Ono et al.||2012 Schenker et al. 


2012 


1. 



The redshift evolution in the EW distribution at z > 6 
is likely to be connected with the reionization process rather 



than with secular evolution process in the galaxies ( Forero- 
Romero et al.|2012| also see |Shimizu et al.| ( |2l)TTj )). Although 
current observations do not allow for robust conclusions, an- 
alyzing EW distributions provides a powerful technique to 
probe the population of star forming galaxies - as well as the 
ionization state of the IGM - during and after reionization 
has been completed (e.g. Dijkstra et al.|20lT also see Dayal 
fc Ferrara|([20T2l)). 



Recently, Fumagalli et al. (20111 showed that at low 
levels of star formation (SFR < 1M@ yr~ ), the effects of 
the statistical description of the stellar Initial Mass Func- 
tion (IMF) and the Cluster Mass Function (CMF) emerge. 
Among the expected effects - gathered under the name of 
'stochasticity' - the most relevant for our discussion is that 
the flux ratio between ionizing and Far UV photons fluctu- 



ates around the mean expected value. Fumagalli et al. ( 2011 1 
show how the clustering of star formation together with the 
sampling of the Initial Mass Function (IMF) cause the ratio 
of the Ha flux and FUV flux density to fluctuate around the 
mean significantly for local dwarf galaxies. This effect should 
be also noticeable in the ratio of Lya flux to FUV flux den- 
sity i.e. the EW. In other words, stochasticity may affect the 
intrinsic Lya EW-PDFs of faint LBGs (-18 <;Muv ^-13). 
Because Lya EW plays an important role in understanding 
high-redshift galaxy populations and constraining the reion- 
ization process, it is important to understand quantitatively 
the impact of stochasticity. 

Furthermore, observational constraints on the 'escape 
fraction' of ionizing photons from star forming galaxies are 
derived from the observed ratio of LyC flux density to FUV 
flux density, which is subject to the clustering of star for- 
mation and stochastic sampling of the IMF, especially at 
SFR <T Mq yr~ . It has long been realized that such galax- 
ies can provide the bulk of the ionizing photons needed to 
reionize the IGM and to maintain it ionized (at a fixed es- 
cape fraction, see e.g. Fig 31 of |Barkana fc LoebpOOlj ) . A 



more recent analysis by Kuhlen & Faucher-Giguere (20121 



echoes these findings. Kuhlen & Faucher-Giguere (20121 



show that in plausible scenarios for reionization, star form- 
ing galaxies with very faint magnitudes M\ im ~ — 13 (corre- 
sponding to SFR ~ 10 _2 MQyr _1 ) contribute significantly to 
the overall budget of ionizing photons. While high-redshift 
galaxies with such low star formation rates remain too faint 
to be detected even by future facilities, the observation of 
galaxies with SFRs= 10 -2 — 1 Mq yr _1 at lower redshifts 
z < 4 can help us to understand the physics of such sys- 
tems. For instance, the escape fraction of ionizing photons 
from galaxies with SFR within the upper end of that range 
have already been constrained by observations at z ~ 3.0 
( Iwata et al.|2009 who constrained the escape fraction down 
to SFR~ IMq yr~ ). The possible dependence of these es- 
timates on stochasticity provides an additional motivation 
to study the consequences of the stochastic sampling of the 
IMF. 

In this paper we set out to quantify the effects of 
stochasticity in the EW distribution of dwarf Lya emitting 
galaxies. We implement highly simplified models for dwarf 
galaxies at high redshift to isolate the stochastic effect. The 
SFH is constant, single metallicity, without any extinction 
effects. We stress that our objective is not to explain or in- 
terpret current observations of LAEs and LBGs, but instead 
simply to isolate and highlight the impact of stochasticity 
on predicted values for the Lya EW. 

This paper is structured as follows. In !|2]we summarize 
the expected effects on the equivalent width stochasticity 
and present the numerical experiments we perform to quan- 
tify this effect. In [J3]we present the raw results that help us 
in Sj4]to asses the impact of stochasticity in different aspects 
of the interpretation of observations of LAEs and LBGs at 
high redshift. Finally, we summarize our conclusions in ij5] 

Throughout this paper we will study the 'Lyman con- 
tinuum' (A ~ 700 A), ionizing continuum (A < 912A), and 
Far UV (FUV) continuum (~ 1500 A), where all wavelengths 
are rest frame. Whenever we refer to UV we imply FUV. 
Also, the Lya equivalent widths are measured in the galax- 
ies' restframe. 



2 STOCHASTIC EFFECTS IN THE 
EQUIVALENT WIDTH 

To be consistent with observations, we define the Lya rest- 
frame equivalent EW as the ratio between the intensity of 
the Lya line to the flux density redwards of the lineQ 



EW 



L 



Lya 



At 



-'Lya 



L\,uv ^uv Lv,uv 



(i) 



1 Formally, the equivalent width is measured relative to the flux 
density in the continuum just redward of the Lya resonance, 
i.e. at A ~ 1216 A+e. The flux density at this wavelength is 
~ (Ac/Atjv)' 3 times larger than at Ajjv = 1540 A, where /3 ~ — 2 
punlop et al.|2012] |Finkelstein et al.|2011[ | denotes the observed 
slope of the UV continuum at the redshift and UV-magnitudes of 
interest. Therefore, if we had corrected for the slope in the spec- 
trum of t he UV continuum , then we would have found EWp ~ 70 
A for the|Kennicutt| |l998 ) star formation calibrators (e.g. |Dijk-| 
stra fc Westra|2010| . 
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where subscript v or A indicates the flux density in frequency 
or wavelength, respectively. The wavelengths and frequen- 
cies associated with the Lya transition and the UV con- 
tinuum (from the central wavelength for the filter used in 
SLUG) are \ a = 1216A, A uv = 1540A, u Lya = 2.47 x 
10 15 Hz and wv = 1.97 x 10 15 Hz. 

We calculate the intrinsic luminosity of Lya assuming 
that a fraction 0.67 of the ionizing photons are converted 
into Lya. We obtain LLyn=co x Qh where Co = 1.04 x 10 -11 
erg photon -1 ( Schaerer||2003 \ . We then have 



EW : 



Auv coQh 



(2) 



In order to compare results on the same grounds re- 
sults from different metallicities or star formation rates we 
introduce a dimensionless constant 



M 



EW 
EWn 



(3) 



where EWo is the equivalent width expected value for high 
values of the star formation rate when the stochasticity ef- 
fects are negligible. 

We point out that under this convention EWo = 
96 A when considering a 'standard' conversion between 



Lya luminosity, UV luminosity density (Kennicutt 19981, 
which are L^uv =SFRx8.0 x 10 27 erg s" 1 Hz 1 and 
L-Ly a =SFRxl.O x 10 42 erg s _1 . We emphasize that an im- 
plicit assumption in the conversion is a Salpeter IMF in the 
mass range 0.1-100 MQfor solar abundance. In our numeri- 
cal experiments, EWo = 110 ± 12 A, which is a result of the 
lower metallicities and the higher upper limit in the IMF 



mass range assumed in SLUG (details are given in 1 2.2 J. 
These minor systematic corrections due to different condi- 
tions in stellar metallicity, IMF or star formation history 
won't affect our main conclusions. Throughout this paper, 
we will adopt EWo = 110 A. 



2.1 Stochasticity 

The description of star formation is probabilistic. At its core 
we find the concept of the Initial Mass Function (IMF), 
4>{m), which describes the relative abundance of stars of a 
given mass to be formed. The IMF can be mostly determined 
by two pieces of information. The first one is the mass inter- 
val, m m i n — m max that the newly formed stars can have. The 
second is the change of relative abundance for different stel- 
lar masses, which is usually parametrized as <j>(m) oc m 1 , 
if 7 > less massive stars are more probable to form that 
high mass stars. 

For large masses of newly formed stars it is possible to 
guarantee a good sampling of the IMF. For a low rate of 
star formation the poor sampling of the IMF makes possible 
that the rarer most massive stars near m max don't appear 
in a burst of star formation. 

Star formation is also thought to proceed in clusters. 
This provides another important astrophysical element in 
the statistical description of star formation that can be im- 
pacted by stochasticity effects. In analogy to the stellar mass 
function, one can also parametrize the abundance of stellar 
clusters by the cluster mass function (CMF). Low values of 



the star formation rate can also affect the sampling of the 
CMF. 

Furthermore, stars of different mass have different evo- 
lutionary tracks. The inhomogeneous IMF sampling will pro- 
duce fluctuating stellar populations with different luminos- 
ity histories, this can impact the ratio of ionizing to far UV 
photons produced as a function of time. 

The effect of stochasticity has been explored before in 
the production of Ha in dwarf galaxies in the local universe 
( Fumagalli et al.|201l}|Haas fc Anders|20l6"| ), but not on the 
effect on the EW of Lymana emitting galaxies. Our aim is 
to quantify this impact and estimate how current and future 
observational campaigns can be affected by these effects. 

2.2 Grids of SLUG models 

We used SLUG ('Stochastically Light Up Galaxies') a fully 



stochastic code for synthetic photometry of galaxies ( Fuma- 
galli et al.|2011 da Silva et al.|20T2 1 . The code samples the 



CMF and IMF of a star forming galaxy by forming indi- 
vidual star clusters and following their spectral evolution. 
At low values of a time averaged continuous SFR the series 
of bursts associated with the cluster formation start to be 
visible. Furthermore, low SFR values imply a low number 
of stars to sample the IMF, lowering the probability of find- 
ing massive stars. All the effects of finite sampling in mass 
and time have an impact on the instantaneous strength of 
the Lymana emission line. SLUG assumes that a fraction f c 
of the total stellar population is formed in clusters. In that 
case, the clusters also sample a cluster mass function(CMF) 



V(M e , 



AC 



The stars in the cluster are then populated 



according to the selected IMF. 

In this stochastic framework, a galaxy with a time av- 
eraged constant SFR does not level to a constant SED, but 
instead it continuously fluctuates in time. For a sample of 
independent galaxies with the same average SFR one could 
thus associate a distribution of line and continuum intensi- 
ties. 

All models used in this paper are based on a Salpeter 
IMF in the mass intervals 0.08-120 M Q . We consider two 
kinds of models of varying metallicity Z = 4.0 x 10 -4 , 4.0 x 
10 -3 (which corresponds to Z = 0.02Z Q , Z = 0.2Z Q ). We 
assume that all the stars are formed in clusters, f c = 1. 
Assuming unclustered star formation, f c — 0, reduces the 
scatter in the flux fluctuations we report here by about an 
order of magnitude ( |da Silva et al.||2012[ ). The CMF is de- 
scribed by the cluster mass interval is 20 — 10 7 Mq and the 
slope is ft = 2.0. The stellar libraries are Padova asymptotic 
giant branch (AGB) tracks. The stellar Spectral Energy Dis- 



tributions (SED) correspond to Lejeune et al. ( 1997 1998 1 



where the Wolf-Rayet stars treatment follows |Smith et al.| 
p002| ) and |Hillier fc Millerj ( fl998t > . We run 500 galaxies for 
each metallicity and star formation rate comprised in the 
range 10~ 3 -1 Mq/jt. Table 1 lists the exact values for the 
SFR and metallicity that were used. Each model runs for a 
period of 200 Myr with outputs every 1 Myr. Running for a 
period longer than 100 Myr ensures a stable base level for 
the mean FUV and ionizing fluxes. 

We use the last 10 outputs per galaxy to slightly in- 
crease the statistical sample. Being pessimistic, one can as- 
sume that each simulated system contributes only 1 inde- 
pendent datum instead of 10, from that point of view, one 
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Figure 1. Different Probabiiity Density Distributions for M., the 
normalized EW. Each symbol corresponds to the binned SLUG 
results for a different SFR=10- 3 , 1CT 2 , lO" 1 and lMQyr" 1 and 
a metalliciticy f 0.20Zq. The lines are the best fit of Eq. Q to 
the points. 
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Figure 2. Fraction of galaxies as a function of SFR that have an 
EW larger than a fixed value EW C . Each line represents different 
threshold in EW C = 20, 30, 40 and 50 A . This figure quantifies 
the impact of the stochasticity on the fraction of confirmed LAEs 
based on a minimal EW detection threshold. 



can be sure that all the results for the PDF and derived 
quantities are accurate within 1/500 ~ 0.02. In other words, 
one might be concerned that the 10 outputs of the same 
galaxy may not be independent and might bias the PDF. To 
quantify this we perform the following test: for 4000 galaxies 
with a star formation rate 8.0 x 10 -3 M0yr _1 we calculate 
the PDF in two different ways. First by taking 10 outputs of 
400 galaxies and second by using only 1 output for each one 
of the 4000 galaxies. We find that the absolute difference 
between the PDF(TM) constructed from the two different 
samples is always on the order of ~ 0.02 and conclude that 
a possible bias is negligible. We perform an additional test 
by taking 1 output and 2 outputs for 4000 galaxies, with 
spectral outputs spaced by lOOMyr during 500Myr and find 
that the difference between the PDF(.M) constructed from 
these samples is always on the order of ~ 0.05 and con- 
clude that a possible bias from the chosen timesteps is also 
negligible. 

For each galaxy we obtain the flux per unit frequency 
i^FUV at a wavelength of 1540 A, together with the ionizing 
photon rate Qh- These quantities are used to calculate the 
EW as described by Eq |2| . The value for EWo is calculated 
as the arithmetic average EW for 30 models of galaxies with 
a SFR of 5 M Q yr _1 , we obtain EW = 110 ± 12 where the 
error quoted is estimated from the standard deviation from 
the mean. 



3 RESULTS 

The net effects of stochasticity introduce a SFR dependent 
scatter around a central value. We describe this effect as a 
Probability Density Distribution (PDF) for the M variable. 
We parametrize the PDF by a double power law in order 
to better quantify how it depends from changes in the SFR 
and the metallicity, 



Figure [T] presents the results for the PDF P(jW|SFR), 
the symbols represent the results from the simulations and 
the lines the fit to those points. Given the limited number 
of models the PDF is reliable only for values P(.M|SFR) > 
10 -2 , likewise the deviation of the fit with respect to the 
simulation points is on the order of 1 x 10 -2 . The star for- 
mation rate enters through the constants Mo, a, 7, and Pq. 
The expected result for a perfectly sampled IMF is a delta 
function around M = 1 (Leitherer et al. 19991, while our 



P(M\SFR) = P 



M 

Mo 



+ 



M 

Mo 



(4) 



results show a wide distribution. Each line in Figure [T] cor- 
responds to a different SFR= 1 x 10~ 3 , 10~ 2 , 10" 1 and 1.0 
in units of Moyr -1 for a metallicity of 0.02z7q. The best fit 
values together with their uncertainty from the fitting proce- 
dure for all the SLUG models are listed in Table 1. In these 
values we get a precise quantification on how the PDF shape 
changes with the SFR. From these results we confirm that 
the effect of the metallicity is a second order effect shaping 
the EW PDF when compared to the influence of the SFR. 

For high values of the SFR> 10" 1 M yr _1 the low end, 
M < 1 of P(M) increases sharply compared to the slower 
decrease in the high end M > 1. In this regime the peak of 
the distribution is around the expected value of M = 1 for 
a perfectly sampled IMF. 

This quantitative behaviour presents consequences for 
the interpretation of high EW systems. On the upper end, 
on can achieve values close to what is expected to be limit of 
star forming galaxies: 240 A which corresponds to M = 2.2 
in our models. For SFR< 10 -2 M0yr _1 , the low end of the 
PDF shows a bulkier shape that confirms the existence of 
systems with M < 1 purely as a consequence of their intrin- 
sic star formation stochasticity. We quantify and elaborate 
in the next section the implications of these results. 
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Table 1. Best fit values for the parameters used in the conditional PDF P(A / 1|SFR) (see Eq[4jl. 



SFR (Mgyr" 1 ) 


Z/Z Q 


Mo 


a 


7 


Po 


0.0010 


0.02 


0.75±0.03 


0.13±0.03 


2.67±0.10 


1.15±0.05 


0.0015 


0.02 


0.71±0.02 


0.44±0.03 


2.73±0.09 


1.44±0.04 


0.0025 


0.02 


0.71±0.03 


0.74±0.06 


2.82±0.16 


1.66±0.07 


0.0040 


0.02 


0.76±0.01 


1.03±0.03 


3.49±0.09 


1.86±0.03 


0.0063 


0.02 


0.80±0.01 


1.21±0.02 


3.83±0.05 


1.97±0.01 


0.010 


0.02 


0.77±0.01 


1.81±0.10 


4.01±0.10 


2.37±0.06 


0.015 


0.02 


0.75±0.01 


2.15±0.08 


3.95±0.11 


2.57±0.03 


0.025 


0.02 


0.75±0.01 


2.79±0.08 


4.12±0.09 


2.89±0.03 


0.040 


0.02 


0.79±0.01 


3.18±0.06 


4.66±0.06 


3.09±0.02 


0.063 


0.02 


0.76±0.01 


4.31±0.07 


4.52±0.06 


3.43±0.01 


0.1 


0.02 


0.77±0.01 


5.12±0.13 


4.57±0.10 


3.64±0.03 


1.0 


0.2 


0.85±0.01 


12.49±0.36 


6.01±0.11 


5.23±0.05 


0.0010 


0.2 


0.64±0.05 


0.013±0.040 


2.25±0.15 


1.20±0.10 


0.0015 


0.2 


0.58±0.04 


0.39±0.05 


2.38±0.14 


1.62±0.09 


0.0025 


0.2 


0.64±0.02 


0.62±0.04 


2.73±0.12 


1.72±0.06 


0.0040 


0.2 


0.78±0.02 


0.70±0.05 


3.60±0.18 


1.74±0.06 


0.0063 


0.2 


0.71±0.01 


0.70±0.05 


3.62±0.14 


2.10±0.05 


0.010 


0.2 


0.69±0.01 


1.52±0.07 


3.96±0.15 


2.50±0.05 


0.015 


0.2 


0.72±0.01 


1.85±0.09 


4.34±0.18 


2.61±0.06 


0.025 


0.2 


0.64±0.01 


3.08±0.11 


3.90±0.09 


3.30±0.03 


0.040 


0.2 


0.70±0.01 


2.87±0.13 


4.50±0.17 


3.26±0.05 


0.063 


0.2 


0.71±0.01 


3.46±0.16 


4.87±0.19 


3.58±0.05 


0.1 


0.2 


0.70±0.01 


4.42±0.13 


4.69±0.10 


3.93±0.03 


1.0 


0.2 


0.76±0.01 


11.52±0.43 


6.08±0.15 


5.87±0.07 



4 DISCUSSION 

4.1 High EW outliers and apparent Lymana 
escape fraction 

The ensuing discussion will focus on the interpretation of 
high and low EW systems that naturally emerge through 
stochasticity. High EW vales might lead one to propose the 
presence of very low metallicity (or population III) stars 
while low EW values could cause some galaxies not to be 
selected as a LAE in very deep narrowband surveys. 

A common assumption in the search for galaxies with 
Population III stars is the following: there is an upper EW 
value to consider a star forming galaxy to be normal. This 
threshold value is taken to be EWth ~ 240 A but, as we 
show in the previous section, galaxies with similar or larger 
EW values can be expected by stochasticity effects without 
demanding Pop III stars. The threshold value EWth = 240 
A corresponds to M = 2.2 (Eq. [3]). We can then quantify 
the fraction of such galaxies that have EW that exceed the 
threshold by using Eq. Q. For galaxies with star formation 
rates in the range 0.1-1 Moyr -1 we find that there would be 
between 0.5 — 2.0% of systems with EW values larger than 
EWth- We calculate same fraction directly from the 500 x 10 
points used to build the PDF and find that between 1% to 
3% among them have an EW larger than EWth- 

This suggests that a minimal threshold in EW is not 
enough to assert that an individual galaxy presents signa- 
tures of population III stars. However, also the effects of 
stochasticity cannot be observationally evaluated from the 
measurement of an individual galaxy. Only statistical stud- 
ies of galaxy samples are an appropriate tool to find clear 
signposts either of special stellar populations or stochastic- 
ity. 



Narrowband surveys for LAEs often employ color cuts 
that can be approximated by requiring that the EW of a 
galaxy exceeds threshold value, which typically lies in the 
range EWcut = 20 - 65 A (see §[]}. The result of previous 
sections show that low values of the SFR< 10~ 2 Mq yr _1 
the EW can be lower than half of the expected EWo. In 
a survey that adopts color-cuts that correspond to a high 
EWcut, one could simply miss those objects due to their 
low intrinsic EW. This shows how stochasticity can induce 
selection bias in LAEs observations. 



Not surprisingly, the impact of stochasticity is highly 
SFR dependent. The strong dependence of the shape of the 
.M-PDF on SFR causes a larger fraction of systems to be 
found below a fixed A4 threshold for detection as the SFR 
decreases. In Figure [2] we show the fraction of LAEs that 
are missed as a function of the star formation rate, for dif- 
ferent cuts in the EW expressed by EW C . These fractions 
are calculated directly from the SLUG models, providing an 
uncertainty on the order of ~ 0.01 in the reported fraction 
of detected galaxies, the discrepancy of calculating the same 
fractions from the fitting function is in all cases smaller than 
~ 0.02. The value of EW C = 20 A corresponds to a common 
cut used in present day observations (e.g. Ouchi et al.|2008 1 



A cut of EW C = 50 A represents the cut used in the early 
observational attempts to find Lyman a emitting galaxies 
(e.g. |Rhoads et al.p OOO). In § [O we investigate how such 



a bias may affect predicted number densities of LAEs in 
simplified (but often employed) models of LAEs. 
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4.2 The connection between the LAE-LF and the 
UV-LF 

In this section we study the potential impact of stochas- 
ticity on the predicted Lya luminosity functions & equiv- 
alent width distribution of narrowband selected galaxies 
(i.e. LAEs). We compare two very simple models to iso- 
late and highlight the impact of stochasticity. In the first 
model we ignore stochastic ffects. We further assume that 
each star forming galaxy has a Lya emission line with 
EW=EW = 110 A ( see 8 2h. In the second model, we in- 
stead assume that the Lya emission line of a star form- 
ing galaxy is drawn from a PDF which is set entirely by 
stochasticity effects. We therefore assume that the function 
P(EW|SFR) is given by Eq|i]for the second model. Finally, 
we further assume that there is a one-to-one correspondence 
between SFR and Afuv (given by the standard conversion 
factor between SFR and UV continuum flux, see §[2f, and 
thereforsQ that P(EW|SFR) = P(EW|Af UV ). Under these 
assumptions, we can predict the LAE Lya luminosity func- 
tion, *(£«), for both models as ( |Dijkstra fc Wyithe|2012 I 



*(L Q ) = lnlO J dM vv 4>{Mw) X EW C x P(EW C |M UV ), (5) 

where ^(L a )d\ogLa denotes the comoving number density 
of LAEs in the luminosity range logL Q ± dlogL a /2. Fur- 
thermore, (f>(Muv) denotes the comoving number density of 
star forming galaxies with absolute UV-magnitude in the 
range Afuv ± dMuy/2. The equivalent width EW C denotes 
the EW which corresponds to a Lya line of luminosity L a 
for an absolute UV magnitude Muv (see Dijkstra & Wyithe 
2012 for details). We focus our analysis on z = 5.7 and 
compare to the observed LAE Lya luminosity functions at 
this redshift by |Ouchi et al.| (|2008[). W e take the observed 
0(Afuv) at z ~ 6 from Bouwens et al. (2007) and 'rescale' 



it to z — 5.7 assuming that M uv evolves with redshift as 
M uv = -21.02 + 0.36(z - 3.8) ( |Bouwens et al.|2008[ ). 

Figure |4.2| shows the luminosity function for the first 
model (no stochasticity) as the black solid line, and the sec- 
ond model (with stochasticity) as the red dashed line. The 
data points represent observations by Ouchi et al. (20081. 



Both models are remarkably similar. The reason for this sim- 
ilarity is non-trivial: stochasticity decreases the 'most likely' 
EW for star forming galaxies compared to models that do 
not include this effect (see Fig[T]). This decrease alone would 
reduce the predicted Lya luminosity functions. However, it 
is compensated for by the fact the A4-PDF has a tail to 
values A4 > 1. The presence of this tail allows for an ad- 
ditional contribution to $(I Q ) (at a given L a ) by UV-faint 
galaxies that would not be able to contribute in the ab- 
sence of stochasticity. Our prediction further assumes that 
P(EW|SFR > 1M Q yr" 1 )=P(EW|SFR= 1M yr" 1 ), and 
therefore likely overestimates the impact of stochasticity at 
SFR> 1M© yr -1 (which corresponds to L a > 10 42 erg s" 1 ). 
We therefore safely conclude that stochasticity barely im- 



2 Of course, in reality stochasticity introduces a dispersion in the 
value of SFR at a given Muv, an( l the proper EW-PDF for a fixed 
M uv is P(EW|A/ UV ) = fd SFR P(EW|SFR)P(SFR|M UV ). Be- 
cause the PDF P(EW|SFR) varies relatively weakly with SFR, 
this proper calculation barely changes our results at all. 
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Figure 3. This Figure shows predicted Lya luminosity func- 
tion of LAEs for simplified models with (without) the effects of 
stochasticity included as the red dashed (black solid line). Also 
shown are the observations taken from Ouchi et al. (2008, black 
filled circles). The main point of this plot is that stochasticity 
barely affects the predicted luminosity functions. Both models 
predict more LAEs than is observed, as we have not attempted to 
correct for dust extinction and/or intergalactic scattering, which 
would both reduce the predicted luminosity functions. 



pacts predictions of Lya luminosity functions. This result is 
not sensitive to the uncertainties in the fitting functions for 
the PDF(A / (jSFR) used here. Also shown for completeness 
is the observed z = 5.7 LAE luminosity function from Ouchi 
ct al. (2008). Both models clearly predict more LAEs than 
is observed. This is not surprising as we have not attempted 
to correct for dust extinction and/or intergalactic scattering. 
Both these effects likely reduce the predicted Lya luminos- 
ity at a fixed number density of objects, which would bring 
down the predicted luminosity functions. 

The EW-PDF for a sample of Lya-selected galaxies is 
given by (Dijkstra & Wyithe 2012) 



logi Q 



P(EW|logP Q , mm ,logL Q , max ) = 
P(EW|M uv ,c)^(Muv,c)rflog 10 L a 



(6) 



where P a ,min (Pa, max) denotes the minimum (maximum) 
Lya luminosity of the galaxies in the sample. Furthermore, 
A/uv,c denotes the UV magnitude for which a Lya line of lu- 
minosity L a corresponds to an equivalent width EW (again 
see Dijkstra & Wyithe 2012 for details). We focus on the im- 
pact of stochasticity on the predicted PDF for logP a ,min = 
42 and logL ajmax = 43 (blue dotted line (which corre- 
sponds approximately to 1 Mq yr -1 <SFR< 10 Mq yr -1 , 
and roughly to the range of observed Lya luminosities). For 
comparison, we also show a case with logL Qjm i n = 40 and 
logL ajmax = 42 (red dashed line, which corresponds approx- 
imately to 0.01 M @ yr -1 <SFR< 1 Mq yr" 1 ). The choice 
logL a ,min = 40 is a bit arbitrary, but was motivated by 
the notion that galaxies with SFR~ 10~ 2 M yr -1 can con- 
tribute significantly to the overall ionizing photon budget 
during the Epoch of Reionization (see Kuhlen & Faucher- 
Giguere 2012). 

The black solid line in Figure |4.2| shows the predicted 
EW-PDF for the model with no stochastic effects (which is 
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Figure 4. This Figure shows the predicted EW-PDF for a sam- 
ple of Lyo selected galaxies. The solid black line shows the Dirac 
delta function for the model with no stochasticity effects. The red 
dashed line (blue dotted line) shows the predicted EW-PDF for 
models which include stochasticity effects for Lyc* selected galax- 
ies with L a = 10 40 - 10 42 erg s" 1 (L a = 10 42 - 10 43 erg s _1 , 
see text). Also shown for comparison is the EW-PDF observed by 
Ouchi et al. (2008; the red filled circles represent their photomet- 
ric data, while the blue filled circles represent their spectroscopic 
sample). This Figure illustrates that stochasticity may help ex- 
plain the large observed values of EW of LAEs at z=5.7, although 
it fails to fully explain the observed EW-PDF. 



simply a delta- function). The predicted PDF for the mod- 
els that include stochasticity peak at the lower values of 
EW, and have tails that extend to values that greatly ex- 
ceed EW= 100 A. Also shown is the observed EW-PDF for 
a sample of photometrically selected (red filled squares) and 
spectroscopically confirmed (blue filled circles, taken from 
Ouchi et al. 2008) LAEs. Clearly, stochastic effects may help 
to explain large observed values of EW. However, additional 
processes such as dust extinction and/or intergalactic scat- 
tering, which impacts differently the UV continuum and the 
Lya line ( Forero- Romero et al. 20111, are needed to fully 



explain the observed EW-PDF of LAEs at z=5.7. 



4.3 Stochastic effects on the Lyman continuum 

The distribution of EW at a given averaged SFR is the result 
of fluctuations in the fluxes of ionizing and far UV contin- 
uum (non ionizing) photons. The ionizing photon production 
is dominated by short-lived O & B stars, while less massive 
A stars can contribute significantly to the non-ionizing pho- 
ton flux. Because of the different lifetimes of these stars of 
different masses, the ionizing and non ionizing fluxes - and 
therefore their ratio - vary on different timescales and with 
different amplitudes. 

Simultaneous fluctuations in the ionizing and non- 
ionizing fluxes can thus impact observational estimates of 
the ionizing escape fraction, which is usually estimated from 
a relative escape fraction: 



f (/Lyc//l500)obs , \ / 7 N 

/esc.rel = 77 T? s eX P (TIGm), (7) 

IjLyC/ 71500 Jstcl 

where the (/L y c//i50o)stei is the intrinsic flux ratio for 



a stellar population model, (/Lyc//i50o)obs is the observed 
flux ratio and tigm is the LyC optical depth of the IGM at 
the redshift of emission. 

Based on Subaru/Suprime-Cam imaging |Iwata et al.| 
(20091) claimed a detection of 10 LAEs and 3 LBGs at z ~ 3.1 



with a significant emission in LyC out of a sample of 125 
LAEs and 73 LBGs. Recent measurements seem to con- 
firm such measurement for 8 LAEs using spectroscopy with 
VLT/VIMOS and SUBARU/FOCAS. The results indicate 
a LyC excess with respect to non-ionizing UV at ~ 1500A 
expected by standard stellar populations ( Inoueet al.|2011 |. 
However, observation targeting LyC in LBGs at z ~ 1.3 us- 
ing deep Hubble Space Telescope imaging do not find any 
excess ( Siana et al.|2010 i. 

The sources present an excess of a factor of ~ 2 — 3 times 
larger than the expected intrinsic stellar ratio for /Lyc//i500- 
Some theoretical analysis requires a very top heavy IMF to 
explain the observed ratios or fractions of 10% in mass of 
very young < 1 Myr of extremely meta l poor Z < 1.0 x 
1O _4 Z0 or metal free stellar populations ( Inoue|2010 
et al.|2011| | 



Inoue 



The stochasticity effects present an important element 
to be considered in such analysis: The ratio M can be di- 
rectly interpreted as the excess of the /i y c//i500 ratio with 
respect to the constant expected value. The star formation 
rates for the objects with LyC excess are on the order of 
~1M yr" 1 (Muv £ - 19, see Fig 5 of Siana et al. 2010 
and where we ignore dust in the conversion from Muv to 
SFR). In this range 10 ± 1% of the galaxies will show an 
excess 1.5 — 2.0 times the expected stellar fhyC / fisoo ratio 
only due to stochasticity (Figure 1), this is consistent with 
the fraction of galaxies observed with those characteristics. 

The stochasticity might also explain why the LyC ex- 
cess is not observed in the recent z ~ 1.3 Hubble imaging by 



Siana et al. (20101. The LBGs observed at that redshift have 
systematically restframe magnitudes A/uv < — 20, while the 
LyC excess sources are fainter than this limit. This transi- 
tion roughly corresponds to a SFR ~ 5 Mq yr -1 assuming 
a standard conversion factor, while the fainter systems re- 
ported with an excess have Muv ~ — 18 and SFR ~ 1 Mq 
yr - . The dependence of the stochasticity effects on the SFR 
make it more likely to detect the LyC excess in the fainter 
sample. Probably the stronger effect is of statistic origin, 
given that the observations at z ~ 1.3 target 15 LBGs, while 
at z ~ 3.0 a set of 198 galaxies were observed, this provides 
one order of magnitude less objects to sample the P(A4) 
distribution and find galaxies with a LyC excess. 



5 CONCLUSIONS 

In this paper we have investigated the quantitative impact of 
stochastic sampling of the stellar initial mass function (IMF) 
on the Lya, and ionizing photon emissivity of dwarf galaxies. 
We have used the SLUG code to simulate the spectral energy 
distribution for restframe wavelengths A <c2000 A of galaxies 
with time-averaged star formation rates in the range (10 -3 — 
1) M yr" 1 . 

From our predicted Far UV continuum flux density (at 
(A) = 1500 A) and the ionizing flux density (at (A) < 912 
A) we construct probability distribution function (PDF) for 
the variable A4 that represents the ratio of the measured 
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restframe equivalent width (EW) of the Lya line to the ex- 
pected constant value in the absence of stochasticity (EWo), 
i.e. M sEW/EWo. 

We find that the A1-PDF can be represented by a dou- 
ble power law in M, whose parameters are highly dependent 
on the SFR. We emphasize that the results derived from 
these experiments are sensitive to the degree of the clus- 
tering properties of the stars. In the case of f c = where 
the star formation is completely unclustered, the dispersion 



around the mean is reduced by a factor of 10 (da Silva et al 



20121. The results are also sensitive to the star formation 



history assumed: in the case of a bursty star formation his- 
tory, one could expect the scatter around the mean for M 
to be higher than what we derive here. 

Our results show that it is possible for galaxies to have 
both extremely low and high A4 values - especially at in- 
creasingly low SFR - and investigate implications of this 
result. In particular, we show how the existence of galax- 
ies with low M values can induce an observational bias if 
galaxies are selected by cuts in the equivalent width (as in 
narrowband surveys for LAEs). This bias may cause dwarfs 
with lower SFRs < 0.1 M^yr^ 1 more likely to be missed 
in narrowband surveys. On the other hand, the existence of 
galaxies with high Ml implies that a high EW cannot un- 
ambiguously be interpreted as a defining characteristic of 
primordial population III galaxies (Schaerer 2003 Johnson 



et al.|[2009] |Pawlik et aT^|20TT) >. Finally, we use the frame- 



work presented in Dijkstra & Wyithe 2012 to quantify the 
impact of stochasticity on the predicted luminosity func- 
tion and EW-PDF of a sample of Lya emitters (LAEs, i.e. 
narrowband selected galaxies) at a fixed redshift. We find 
that while stochasticity does not appreciably affect the LAE 
abundance modeling, it may help to explain the large EW 
(EW ^100 - 200 A) observed for a fraction of LAEs. 

Finally, we have also shown that stochasticity can af- 
fect inferred constraints on the escape fraction of ionizing 
radiation, most strongly in galaxies that form stars at rates 
SFR <1Mq yr" 1 . Given that such galaxies are thought to 
dominate the reionization process at z > 6, it is important 
to consider the effects of stochasticity when using 'local' 
observations to put constraints on their ionizing luminosi- 
ties and escape fractions. In particular, the observed anoma- 
lies in the ratio of LyC to non-ionizing UV radiation in so- 
called 'Lyman-Bump' galaxies at z — 3.1 reported by |Iwata| 



et al. (2009| can be explained naturally in the context of 
the stochastic sampling of the IMF, without resorting to 
unusual stellar populations. Moreover, the absence of such 
objects in a sample of 15 LBGs at z ~ 1.3 (ISiana et al 



2010) is also expected as these galaxies are forming stars at 



rates at which the stochastic sampling of the IMF becomes 
negligible. 
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